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Abstract 

Abstract: The next-to-leading order (NLO) evolution of the parton distribution 
functions (PDFs) in QCD is a common tool in the lepton-hadron and hadron-hadron 
collider data analysis. The standard NLO DGLAP evolution is formulated for inclu- 
sive (integrated) PDFs and done using inclusive NLO kernels. We report here on the 
ongoing project, called KRKMC, in which NLO DGLAP evolution is performed for 
the exclusive multiparton (fully unintegrated) distributions (ePDFs) with the help 
of the exclusive kernels. These kernels are calculated within the two-parton phase 
space for the non-singlet evolution, using Curci-Furmanski-Petronzio factorization 
scheme. The multiparton distribution, with multiple use of the exclusive NLO ker- 
nels, is implemented in the Monte Carlo program simulating multi-gluon emission 
from single quark emitter. High statistics tests (~ lO"*^*^ events) show that the new 
scheme works perfectly well in practice and, at the inclusive (integrated) level, is 
equivalent with the traditional inclusive NLO DGLAP evolution. Once completed, 
this new technique is aimed as a building block for the new more precise NLO parton 
shower Monte Carlo, for W/Z production at LHC and for ep scattering, as well as a 
starting point for other perturbative QCD based Monte Carlo projects. 
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We report on the ongoing effort on the exclusive Monte Carlo (MC) modeling of 
DGLAPjl] evolution at the NLO level using classic work of Curci Furmanski and Petronzio 
(CFP) [2] as a guide and reference. 

The so-called factorization theorems [3l HI [5] in Quantum Chromodynamics (QCD) 
are stating that in the high energy scattering process of hadrons, with an experimentally 
identifiable single large scale (effective mass, transverse momentum etc.) one may reorga- 
nize the infinite order perturbative expansion in terms of Feynman diagrams, such that 
all coUinear (mass) singularities are encapsulated into certain well defined objects, called 
parton distribution functions (PDFs) or parton fragmentation functions (PFFs or jets), 
while the remaining part, free of such singularities, forms the so-called hard process part 
(coeff. function). The soft singularities due to zero mass gluon emissions are shown not 
to disturb or invalidate this picture [4j, if they are averaged/integrated over the phase 
space and properly combined with the virtual contributions. In the physical gauge the 
PDF/PFF part consists of a well defined Feynman diagrams with the ladder topology. In 
the early stage of formulating practical QCD perturbative methodology it was found that 
the most economical way of dealing with the PDF/PFF parts of the process was to define 
them as inclusive as possible, integrating over transverse momenta and summing up over 
all partons emitted from the ladder, keeping control only on the total energy (light-cone 
variable) of the parton entering the hard process, and its type. Such inclusive (collinear) 
PDF is widely used until today in most of practical QCD calculations, especially for the 
initial hadrons. 

The only exception is the so-called parton shower Monte Carlo (PSMC), where one 
gains access to all momenta and other quantum numbers in PDF/PFF (ladder) part of 
the hadronic scattering, for every incoming hadron or outgoing jet. Originally the main 
role of PSMCs was to describe hadronization of the partons, but they have gradually 
absorbed the leading order perturbative QCD (pQCD) description ^ ^Tj of the ladders 
(PDFs, PFFs). With the growing sophistication of the high energy (HE) experimental 
detectors PSMC became indispensable for understanding data in any modern experiment. 
However, although pQCD calculations using inclusive PDFs have evolved enormously in 
their sophistication (evolution of PDFs at NLL, NNLL level, corrections to hard process 
at NLO, NNLO, new ingenious methods of calculating tree-level multiple parton distri- 
butions and more) the PSMCs have stayed, from the pQCD point of view, where they 
were 25 years ago, that is at the (improved) LO/LL level, until today. This lack of the 
progress is not fully understood, but most likely the main reason was that computers fast 
enough were not available and due to the difficulties in reformulating QCD factorization 
theorems into a form suitable for stochastic simulation (MC) methods. 

We are reporting on the first serious attempt to upgrade parton shower MC for a single 
incoming quark (non-singlet PDF) to the level of the complete NLO|^ This will be done, 
as in early days of pQCD, in the physical gauge, including first order real and virtual 
corrections to the basic ladder describing LO level showering of one incoming quark. 
Our Monte Carlo implementation of NLO DGLAP evolution is: (1) based firmly on 

^See also refs. [Uinilin] for similar effort in this direction. 
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Feynman diagrams and standard LIPS, (2) based rigorously on the collinear factorization 
(eg. EGMPR [3j), (3) implementing exactly NLO MS DGLAP evolution at the inclusive 
level, (5) defining fully unintegrated exclusive ePDFs (the integrand of inclusive PDFs), 
(6) performing NLO evolution by the MC itself (no use of backward evolutionjG]). 
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Figure 1: Raw factorization theorem in the physical gauge. 
Factorization scheme of EGMPR [3], see Fig. [ij was customized to MS by CFP [2]: 



1-Ko V W V 'e 
T(a,-] = (-^] =1 + K + K(^K + K®K®K + ..., (1) 



e 



1-K 



K = PKq ■ , C = Co- ^ 



l_(l_P).Ko' ^ l-(l-P)-i^o' 

where the ladder part F corresponds to MC parton shower and C is the hard process part 
NLO kernels were extracted in ref. ^ from the coefficient of a single pole ^ in T. Projection 
operator of CFP, F 

— Pspin Pkin PPy consists of the kinematic (on-shell) projection operator 
Pkin, spin projection operator Pspin and PP extracting pole part ^, > 0. 

In our MC solutions we use the standard interpretation of the collinear e-poles: 

1 /-^^dF fk^y 







However, the ladder part in CFP/EGMPR scheme features enormous cancellations, as can be 
seen already at the LO leveQ 



1 - ( 1 - e"^ 



--l+[l-e-^) + [l-e-^) +..., (3) 

while from RGE and explicit LO calculation we obtain readily F = e'^^ = 1 + ^ + + ... In 
the MC we need this exponent manifestly, if possible directly from the Feynman diagrams! 

The above exponential nature of the QCD evolution of PDFs is manifest in the following 
master formula 



^ = Co • KlKo] ■ expTO {'^0 • K[Ko]} 



^ = ^° 1 K 



l-Ki 



Multiplication symbol • means full phase space integration (f-k while convolution ® only the inte- 
gration over the 1-dim. lightcone variable. 

"^Omitting for simplicity 1/e" poles due to running of the coupling constant from the consideration. 
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which is serving as a generating functional of the exclusive parton distributions implemented in 
the MC. Here, exp^^Q is the time ordered exponential in the time evolution variable t = In /x, 
where is factorization scale variable. Operator ^ is defined recursivel}j^ 

%^[Ko ■ Ko ■ Ko] =Ko-Ko-Ko- K{''Ko} ■ KA"'Ko} ■ F'.J^^i^o} - • • • 



More terms in the recursion can be obtained by expanding LHS and RHS of eq. ([l]) in powers of Kq 
and using definition of exp^Q (see below). The key point is the definition of modified projection 
operator (a) it does spin projection as F of CFP, (b) it sets its incoming momentum on-shell 
in the part of the diagram towards the hard process, (c) it acts on the integrand of the Lorentz 
invariant phase space (LIPS), before integration, (d) it sets upper limit on the phase space 
for all its own real (cut) partons, eg. fj, > s{ki, ..,kn) = max(A;J'), (e) our preferred choice is 
rapidity ordering choice; s{ki,..,kn) = a{ki, .., kn) = ma-x.{kf / ai) , at = kf/E, (f) ^'^{A) acts 
on A which is at most single-log (collinear) divergent and extracts this singularity from the LIPS 
integrancj^ (g) ^'^{Kq) is legal, as Kq is single-log divergent, (h) nesting like ^ ' [Kq - {l — ^ ' (Kq))] 

is allowed, as long as its argument is no more than single-log divergent, (e) does not include 
PP operation. Finally, the time ordered exponential reads: 

expTO (f;{^})^ = 1 + F;{A} + f;{«M} • K^rA} + F'^i^^A} • K^^A} • K^^A} + . . . (6) 

where notation {^A} means that s = a(ai,...,a„) = max(ai, a„). For instance for n = 3 
the entire integrand multiplied by ^^>s3>s2>si- Variable /x is constant, while Sj depend on the 
4-momenta integration variables. 

Master formula of eq. is very important for us, as it serves as a generating functional 
of the exclusive distributions beyond LO implemented in the MC. In CFP scheme the time 
ordered exponential is present for the inclusive PDFs, and results from the renormalization group 
equations, not directly from Feynman diagrams. Deriving eq. (Q directly from diagrams at any 
perturbative order remains an open important problerrj^ For the time being we check order 
by order that: (i) at the inclusive level we maintain full compability with the CFP scheme, (ii) 
we reproduce fixed order matrix squared element times LIPS results wherever possible and/or 
necessary. 

From now on, in the factorization formula (|4|, we focus on the exclusive PDF (ePDF) 
T>(fi, ki, kn), which is the integrand of the inclusive PDF: 

D{f,) = expTO {P {'Ko ■ K[Ko]}) {ii) = e^VTo{K), (7) 

The x-dependent PDFs (inclusive) are obtained by means of inserting 5{x — x{ki, k^)) in the 
integrand]^ D{fj,) — > D{fi,x) = D{^)x. 



^ Similarly as /3- functions in Yennie-Frautschi-SuurafT? subtraction scheme. See also |13| . 
^For instance by rescaling all kj — ^ XkJ and taking coefficient in front of 1/A term. 
^In particular the effect of the running coupling costant deserves careful discussion. 
^Such an insertion we shall often mark as {■■■)x- In the MC it means histogramming of x. It also helps 
eliminating x — from the consideration. 
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The standard inclusive PDF, D(fj,, x) = J dLips "Ddi, ki, kn, x), is obtained form ePDF by 
the phase space integration. It obeys by construction the ordinary evolution equation 



df,D{fi,x) = y(^D{fi,x) 



(8) 



with the inclusive DGLAP kernel 
d 



7{x) 



dln{fi) 



{K^)x = j rfLips 5 



-t; ['K„ . 6.[/foi} . (9) 



^ J (iLips 

The LO and NLO truncations of the evolution kernel are: Kj^^ = t'^ {'^Kq}, taken at 0{a^) 
and K^^^ = ^'^ ^'[Kq + Kq ■ (1 - • i^o)}, truncated at 0{a^). The 2PI kernel Kq of CFP 
scheme (non-singlet bremsstrahlung) at LO+NLO is: 



(10) 



where dashed lines are gluons, blobs marked "Virt" may include several (one loop) subgraphs. 
First two terms in the x-dependent T.O. exponential with LO kernel read 



J \ L 

I J" j 



5x=\ + In 



Q 2CFas ( l + x 



qo TT V2(l-x) 



Sx=l + In — Tqqix), 

qo 



where go is IR cut-off and Tgg(x) = ^^f^ (^i^; ' 



The LO exclusive distributioij^p^i 



2 ^^Gai>s resides inside 



J i L 

: g" 1 



(fki 
'214 



>ai>qo PBlrih) = / aidaid a.id(t)i Oq^a^^qo PBlrih 



(12) 



where = ki/ai and ai = |ai| ~ polar angle of the gluon with transverse momentum ki. The 
trivial phase space integration gives Sudakov double log or LO kernel: 



In — I In - 



6 4 



J i L 



In ^ y,,(x)0i_,<,. (13) 
Qo 



In the TO exponential D{Q) = exp^-Q (Fq{JC^'^}) the virtual contribution ^~p^= —Sisr 
factorizes off: 

D{Q) = 1 + F'q{Ko} + F'grKo} ■ KA^'Ko} + P'qV'Ko} ■ K^V'^o} ■ KV'Ko} + ... 



exp I K ^ { ^ 



I i I _l \ L, 





...if 







(14) 



e-5:sa|l +^]^^y _^ n(fe,) pi^?(^,)0a,>a,_,| = e'^ISH {e+^ISR.} = 1. 



'Here 6 is another IR cut-off parameter, the same as in CFP. 
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The x-dependent version of the same D{Q,x) = exprpQ ^Fq{X^'-^}^ reads 



oo / n 



^ da, 



„=i Vj=i -^^-i "J -^0 



1-5 



dZj 



2^ difj 



2^ ^gg(^i) )'^^=a^» 



(15) 



From now on we enter NLO world. The NLO kernelJC^^*^ = ^'^^''{Ko + Kq ■ {1 - ^') ■ Kq) } 
with Ko of eq. (10) (nonsinglet bremsstrahlung) is inserted into T.O. exponent of the NLO ePFD: 



Z^W(Q) = exp(-5Wj (l + F'q{K^} + F'q^-K^} • K,rK^}+ 



(16) 




Where the zero-real-emission part of the kernel — -S"!^^ ^ (wave function renormalization 
up to second order) factorizes and exponentiates. 



The remaining part of the NLO kernel = (| | j — J+j.^.^.^ consists of the 



one-real-emission part including NLO virtual one loop corrections K^"" = ^■ ■{^ — (^TTpT^ 



2K 



and the pure two-real emission part K 



2r 



m 




Eq. (16) 



could be used directly in the MC, if both K^^ and K^^' were positive. Since K^^ is non-positive 
we have to recombine K^^ with 2 real emission distribution, that is to put it back where it came 
from (reversing what the factorization did): 



Z)W(g) = e 



1- 



mm 



:i7) 



m 



In eq. (17) the first row is positive and large (LO + virtual NLO), hence it will go to basic MC, 
while the second row being small and nonpositive can be absorbed into MC correcting weight. 
However, it is not so easy to recombine K^'' with 2 real emissions. The target distribution 

^ P?fy cannot be implemented in the ( 17 ) by means of MC-reweighting ^^j^^ ' , simply 
because it is zero outside simplex Q > a2 > ai > qo, while the target distribution is nonzero in 
the bigger rectangle Q > max(a2,ai) > qq. Going back to original Feynman diagrams we see 
that the above problem turns out to be fictitious, if we properly keep track of the Bose-Einstein 
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(BE) symmetrization: 



if... 


o 


if... 





d^k2 f (fki 



(18) 



2^2 



2)^0 '^Q>max 



(a2,ai)>0 Oa2>aiil3lBik2,ki) + /^IB (A;i , ^2) ) , 



where for purely technical reasons we include an internal ordering Oa2>ai for the already sym- 
metric integrand. See ref. for definition of the NLO 2-gluon function /3ib(A;i, /C2). We call 




Figure 2: Two-real-gluon NLO correlation function. 

sometimes the resulting MC weight the NLO short range correlation function, because it con- 
tributes significantly only if both gluons are non-soft and have transverse momenta (or rapidities) 
almost equal, as seen in the plot of Fig. [2] see |14) . BE symmetrization requires clever reorgani- 
zation of combinatorics, if want to gain on the computation speed by means of excluding terms 
equal zero from the BE symmetrization sum. For instance, BE symmetrization over 3! permu- 
tations of 3 arguments of a single NLO correlation function and one LO spectator distribution 
reduces to only 2 terms 
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E 

{tt} 




(19) 



Permutation tt^ performs the right ordering a(7rj,)3 > '^(77^)2 > ^{Trk)3 ^ given phase space point 
k = {ki, k2, ks). This is denoted in the above by {tt*} = {(123), (213)}, where (123) is identity 
and (213) interchanges (vr^)! and (vrfc)2. Generalization to double insertion of NLO correlation 
function in the ladder in a sketchy graphical form is: 



E E E 

Ja,Jb {K} H} 





The corresponding algebraic formula of the above NLO ePDF can be found in ref. [T 
with the numerical result of the precision MC test. 



(20) 



together 
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For simplicity we are drawing here only of half of the ladder. 
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Figure 3: Feynman diagrams contributing to NLO non-singlet evolution. 



So far we have considered the ~ Cp diagrams, in Fig. [3] The remaining ~ CfCa diagrams 
add new problem: strong cancellations between real and virtual contributions in the NLO cor- 
rectiorp^ due to final state radiation (FSR) Sudakov double log, This enforces exponentiation 
of the FSR already in the LO basic MC, if we aim at positive weight MC events. In the basic 



MC each gluon in the LO ladder is replaced by resolved multigluon 2^ '~'''"^\'-?".'v 



, obtained by 



iterating/exponentiating "soft counterterm" 



in Rf . BE symmetrization is also done for 



FSR gluons, such that the complete NLO correcting weight includes the following sum: 





(21) 



where K-^ 




iQfol 



is the 



CfCa component in the above NLO MC 

weight. Th leading FSR Sudakov double log part is subtracted in X-^. The above sum in the 
MC weight was already tested for single NLO insertion, giving rise to a well behaved weight 
distribution. More testing is under way. 

Summarizing, we report here on a small but essential part of the larger project, with the 
aim of implementing NLO DGLAP evolution of the parton distributions in fully unintegrated 
inclusive form in the Monte Carlo with positive weights (weights equal one). The main result of 
the ongoing study, so far limited to nonsinglet ePDF, is that this is feasible. Once completed. 



Soft limit of this object was analyzed in ref. [H], exposing colour coherence effects. 
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this project will lead to new type of the parton shower MC for the initial state in LHC and other 
colliders with hadron beams. 
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